Molecular Phylogenetics and Evolution

Zingiberales includes eight morphologically distinctive families, but its family-level backbone has remained unstable, especially around Musaceae, Heliconiaceae, Lowiaceae, and Strelitziaceae. We analysed 1566 low-copy nuclear genes from 52 samples, representing all eight families and Pontederia crassipes as outgroup. Concatenated maximum likelihood and multispecies coalescent analyses recovered the same backbone: ((Zingiberaceae, Costaceae), (Cannaceae, Marantaceae)) is sister to (Musaceae, (Heliconiaceae, (Lowiaceae, Strelitziaceae))). Penalized-likelihood dating placed the sampled crown group in the Late Cretaceous, with several deep family-level divergences occurring on short internodes. Analysis of 1248 rerooted gene trees showed that conflict is concentrated on these deep branches and in several shallow clades. HyDe tests of empirical and simulated matrices, each including 62,475 triples, did not support widespread ancient hybridization among the major family-level lineages after filtering against the simulated null model. The nuclear data recover a stable Zingiberales backbone, and the long-standing instability of several deep nodes is best explained by rapid early divergence and extensive incomplete lineage sorting.

The global study of freshwater clams in the genus Corbicula is frequently confounded by invasive androgenetic lineages that experience mitochondrial DNA capture and clonal propagation. In contrast, the endemic Corbicula of Sulawesi's ancient lakes reproduce sexually, offering a uniquely reliable system for mitochondrial population genetics. This study provides the first population-genetic framework for two endemic species, Corbicula possoensis (Lake Poso) and C. linduensis (Lake Lindu), using the cytochrome c oxidase subunit I (COI) marker. We analysed 90 newly generated COI sequences from C. possoensis (six stations) and C. linduensis (three stations), integrated with reference sequences from GenBank, to assess genetic diversity, population structure, and phylogeographic patterns. Hierarchical AMOVA revealed deep divergence between the two lakes ({Phi}_CT = 0.607), consistent with prolonged independent isolation rather than a single shared vicariance event, as the two species do not form a sister pair in the phylogeny. Within Lake Poso, C. possoensis exhibited exceptionally high genetic diversity (24 haplotypes; h = 0.876; {pi} = 0.016) and pronounced micro-geographic structuring into three phylogeographic zones (North: Tentena and Siuri; East: Tando Nceppo and Busogo Beach; Southwest: Bancea and Pendolo), each characterised by distinct haplogroups. Remarkably, the maximum divergence between zones (K2P = 2.33%) approached the interspecific distance between C. possoensis and C. linduensis (K2P = 2.42%), indicating that within-lake mitochondrial divergence has reached near-interspecific levels. Conversely, C. linduensis displayed near-panmixia and extreme genetic depauperation (3 haplotypes; h = 0.246; {pi} = 0.0004), indicating long-term demographic stasis within a restricted habitat. The deep phylogeographic zonation in C. possoensis suggests that its discrete populations should be treated as separate Management Units (MUs) in conservation planning to preserve locally adapted gene complexes, whereas the severely depauperate gene pool of C. linduensis renders it critically vulnerable to environmental disturbance and invasive species, warranting urgent IUCN Red List assessment. To validate these mitochondrial boundaries and inform future conservation strategies, multi-marker and genome-wide reassessments are strongly recommended.

The tribe Spermacoceae (Rubiaceae) comprises a morphologically diverse assemblage of approximately 1,400 species distributed across the Neotropics, Africa, Asia, Australia, and Pacific region. It remains one of the most taxonomically intractable groups in the family, with generic limits repeatedly redefined for more than two centuries. Previous phylogenetic studies based on a limited number of plastid and nuclear markers left numerous relationships unresolved and provided sparse representation of Neotropical lineages. Here, we present the first phylogenomic study of the tribe based on plastome-scale data and expanded sampling of Neotropical taxa. We sampled 121 species representing 55 genera spanning all major clades and generated 123 new plastomes, including 25 species incorporated into a molecular phylogenetic framework for the first time. Maximum-likelihood and Bayesian analyses recovered a highly resolved and strongly supported phylogeny, with uncertainty restricted to a small number of deep backbone nodes. Pollen and seed micromorphology provided additional evidence for evaluating phylogenetic relationships. The resulting phylogenetic framework clarifies generic boundaries across several problematic lineages and supports multiple taxonomic changes. Pervasive homoplasy in seed and floral characters rendered several traditionally recognized genera non-monophyletic, warranting new combinations, including Edrastima oxycoccoides, Stenotis alexanderae, and S. prostrata, and a reassessment of taxa such as Terrellianthus serpyllaceus and Oldenlandia dusenii. We further identify genera requiring additional study and provide an updated key to the 82 recognized genera of Spermacoceae. Together, these results provide the most robust phylogenetic framework yet available for the tribe and establish a foundation for future systematic, biogeographic, and evolutionary research.

Reliable species delimitation underpins biodiversity assessment but remains difficult for organisms with plastic morphology and few diagnostic characters. Multispecies coalescent (MSC) methods can delimit species from genomic data, yet they are rarely tested in taxonomically complex, marine invertebrate groups where they are arguably most needed. We used the three Mediterranean species of the genus Tethya, a rare, well-characterised system within the otherwise taxonomically difficult phylum Porifera-distinguished by multiple independent morphological and ecological characters-to evaluate how robust MSC-based delimitation is in such groups. Analysing 64 single-copy nuclear loci in BEAST2 and BPP, we compared constrained, hypothesis-testing approaches (BFD*, BFdriver, A10) with freer, heuristic ones (SPEEDEMON, A11), and examined their sensitivity to data type, clock model, priors, and the species-collapse threshold. All methods recovered the three recognised Mediterranean species, but the resolution of within-lineage structure was method-dependent. The hypothesis-testing approaches consistently supported six lineages, robustly across data types and model assumptions, whereas the heuristic approaches proved less stable. Configurations without a priori species hypotheses often failed to converge or were computationally intractable, a problem compounded by the relaxed clock. In SPEEDEMON the outcome changed with the collapse threshold. Because our system lacks an independent reference point to calibrate this threshold, any delimitation based on it is poorly constrained. We conclude that constrained, hypothesis-testing delimitation is the most robust and reproducible MSC approach, yielding a quantitative, model-based hypothesis that can be weighed against other lines of evidence to inform taxonomic decisions. By clarifying how these methods behave and how their outcomes should be interpreted, our study offers a practical guide for researchers working on comparably complex systems.

Collembola (springtails) are among the most abundant and ecologically important soil arthropods, representing one of the oldest extant terrestrial hexapod lineages, with a fossil record extending to the early Devonian. Despite their relevance, phylogenetic relationships among the four extant orders (Entomobryomorpha, Poduromorpha, Symphypleona, and Neelipleona) have remained unresolved for over two decades. Here, we present the most comprehensive phylogenomic analysis of Collembola to date, comprising 1,127 single-copy orthologues from 145 taxa representing 19 families. To improve orthology inference, we developed a novel HMM-based filtering pipeline that significantly reduced hidden paralogy in BUSCO-derived datasets. Across multiple dataset configurations, gene-jackknife replicates, and various maximum-likelihood analyses, we consistently recovered Poduromorpha as the earliest-diverging lineage. Coalescent-based methods instead highlighted discordant arrangements characterised by extremely short internal branches and low quartet support, a pattern consistent with pervasive incomplete lineage sorting and reticulate evolutionary history. We further dissected the phylogenetic signal by exhaustively evaluating all possible inter-order topological arrangements, both on the full concatenated dataset and gene-by-gene, to identify the most phylogenetically informative loci. These analyses rejected the great majority of previously proposed hypotheses, narrowing support to only two statistically indistinguishable topologies (T11 and T4), with the Poduromorpha-first arrangement consistently favoured across both site-homogeneous and site-heterogeneous substitution models. Finally, with molecular dating, we estimated the origin of crown Collembola in the Early Devonian, with the diversification of the extant orders in the Carboniferous. Several extant genera were estimated to be older than many currently recognized families, highlighting the exceptional evolutionary persistence of springtail lineages and suggesting that lineage longevity should be considered when interpreting higher-level taxonomic diversity.

This study presents a systematic revision of the suborder Astrophorina for the temperate Pacific coast of the United States and Canada. Major findings include a reduction in the number of species previously thought to range into the region from Japan; validation of most Geodia species erected by Lendenfeld (1910), which were later synonymized by de Laubenfels (1932); the formal description of 10 new species (Poecillastra alaskensis sp. nov., Vulcanella explorata sp. nov., Vulcanella rupta sp. nov., Stelletta cardenasi sp. nov., Stelletta nicolenya sp. nov., Stelletta limuwensis sp. nov., Dercitus (Stoeba) giveni sp. nov., Penares anyapax sp. nov., Penares foxi sp. nov., and Thenea diastra sp. nov.); and one new combination, Penares orientalis comb. nov. Extensive SCUBA-based collection efforts yielded new samples for 11 of the 26 species identified in the region, which enabled an integrative taxonomic approach that combined field photography, fresh material for DNA sequencing, and improved characterization of species ranges and morphological variability in previously described taxa. Illumina sequencing generated complete nuclear ribosomal haplotypes for five species, while Sanger sequencing of the 28S and cox1 loci placed 20 of the 26 species within molecular phylogenies. The use of very short "mini-barcode" amplicons also enabled sequence recovery from historic type specimens up to 137 years old. This study additionally reports the discovery of sponge grounds of abundant, large Geodia at diving depths in Southern California. Together, these results substantially advance our understanding of global astrophorid diversity and systematics, and the biogeography of sponge diversity in the Northeast Pacific. Note about species names: this pre-print is not intended to be a publication of the associated species names for the purposes of zoological nomenclature.

Chemosensory systems play a central role in host detection, feeding behavior, and habitat selection in hematophagous insects. Here, we performed a comparative evolutionary analysis of chemosensory gene repertoires across 13 species of the Chagas disease vector genus Rhodnius. While gustatory receptors (GRs), ionotropic receptors (IRs), odorant-binding proteins (OBPs), and chemosensory proteins (CSPs) remained globally conserved, odorant receptors (ORs) displayed extensive lineage-specific expansions, tandem duplications, dynamic transcriptomic regulation, and recurrent signatures of positive selection. Major OR expansions were observed in Rhodnius robustus and Rhodnius colombiensis, suggesting increased sensory diversification in ecologically heterogeneous lineages. In contrast, conserved GR1 expression supports the maintenance of ancestral sugar-detection pathways despite hematophagy lifestyle. We further found no evidence of the canonical insect CO2-associated GRs, suggesting alternative molecular mechanisms for CO2 perception in Triatominae. Several receptors, including Orco, also displayed shifts in selective constraints between sylvatic and domiciliary species, consistent with sensory remodeling associated with adaptation to domestic habitats. Together, our results identify ORs as the most evolutionarily dynamic component of the Rhodnius chemosensory repertoire and highlight contrasting evolutionary trajectories among chemosensory gene families during ecological diversification and vector adaptation.

Sexual conflict has been hypothesized as a driver of speciation, though its effects are likely heterogeneous across phylogenies and between sexes. The semi-aquatic bug, which inhabits water surfaces across diverse aquatic environments, has long served as a model for studying sexual conflict. While previous studies have focused on rapid antagonistic coevolution and the genetic basis of sexually antagonistic traits, the macroevolutionary consequences of asymmetrical sexual conflict--particularly male-dominated grasping traits versus female resistance--remain largely unexplored. Within the subgenus Pseudovelia, males exhibit pronounced phenotypic diversification in grasping structures, whereas females show modest, clade-specific resistance traits, suggesting male-biased asymmetric conflict. This system presents a valuable opportunity to examine how sexual conflict influences diversification and asymmetrical trait evolution across lineages. Using 204 individuals, representing over half of the subgenus's species diversity, we reconstructed a time-calibrated phylogeny, quantified diversification rates, assessed sexual conflict intensity across clades, and analyzed correlations between sexual trait evolution and diversification. Our results reveal extensive phylogenetic conflict, particularly within the East Asian clade, driven by introgression and incomplete lineage sorting (ILS). Furthermore, we observe significant phylogenetic heterogeneity in both phenotypic evolution and diversification rates. Notably, a male "trait package" enhancing grasping ability likely drives rapid diversification in the recently radiated "South China" lineage. In contrast, grasping traits involving abdominal segment VIII are associated with lower conflict intensity, facilitating greater evolutionary flexibility in female resistance and resulting in lineage-specific counter-adaptations. These findings highlight the heterogeneous dynamics of asymmetrical sexual conflict in shaping diversification and speciation.

In search of recurrent patterns of postzygotic hybrid incompatibility, we investigated the Holarctic 'Old World Swallowtail' (Papilio machaon Linnaeus) butterfly complex, whose members can easily be crossed in the laboratory. How many species this model system comprises remains unclear, as taxa with highly distinctive larvae but uncertain status come into contact with authentic P. machaon subspecies in Southern parts of the Palearctic region. By determining mitochondrial and ITS2 haplotypes within and away from contact zones, we found that these neighboring populations do exchange genes, as expected from F1 hybrids being generally fertile in our laboratory crosses. Nevertheless, recurrent instances of dysregulation of diapause were uncovered in hybrid progenies. In keeping with Haldane's Rule, pupae of the heterogametic (female) sex were either unable to enter diapause or, in reciprocal crosses, unable to resume development after having initiated diapause, whereas F1 males experienced normal, photoperiod-regulated diapause, but passed on abnormal diapause regulation to part of their female offspring when backcrossed. Comparing male and female pupal weights in hybrid progenies provides estimates of parental incompatibility that allow to rank taxa and predict quantitatively the outcome of additional crosses, as expected if the same regulatory system were involved. A survey of the entomological literature confirms that diapausing pupae that cannot resume development ('perpetual nymphs') are a recurrent feature of interspecific lepidopteran crosses. Moreover, of the two parent species of a perpetual nymph, the paternal one generally has fewer broods per year. These observations are discussed in the light of models of the speciation process.

Demographic factors are intrinsically crucial to evaluate species' extinction risk. However, measuring them remains difficult and time-consuming and the use of genomic summary statistics has been advocated to assess the conservation status of a species. In the present study, we estimated (i) the census number (Nc), (ii) effective population size (Ne) over three different time periods, recent, historical and ancient, (iii) neutral genetic diversity ({pi}4), and (iv) a measure of the efficacy of purifying selection ({pi}0/{pi}4) for 101 plant species using population genomic sequencing data. Twenty-one species are from the Plant Species with Extremely Small Populations (PSESP) program of SW China. Threatened species exhibited significantly lower Ne, Nc, {pi}4, and weaker purifying selection, but had a higher Ne/Nc ratio than non-threatened ones. Nc was the main determinant in identifying conservation status, and contemporary neutral genetic diversity was predominantly influenced by historical Ne. In the absence of demographic information, genetic parameters are a good proxy of conservation status, likely because currently threatened species also had a low historical population size. In summary, our findings suggest that direct estimates of Nc are more useful than {pi}4, although the latter remains a valuable conservation indicator. Hence, efforts such as the PSESP should be extended.

During mycological surveys in Yunnan Province, China, specimens of a fungus producing massive, upright stromata up to 50 cm high and individually 2.2 Kg in weight were sampled. Through an integrative taxonomic approach combining detailed morphology, multilocus phylogeny (ITS, LSU, RPB2, TUB2), and phylogenomic analyses, this fungus is proposed as the new species Dianjunus rex gen. et sp. nov., the type of the new family Dianjunaceae (Xylariales). Phylogenetic analyses robustly place Dianjunaceae as a distinct sister clade to Graphostromataceae. Divergence time estimation dates the origin of this family to the early Paleocene (~65 Mya), coinciding with the post-K-Pg extinction period, when an estimated 75% of all plant and animal species went extinct, and a significant ecological reorganization of life on earth happened. The stromata of D. rex represent the largest fructifications documented within the Ascomycota, significantly expanding the known morphological range of the Xylariales. The study provides a comprehensive description, including a nodulisporium-like anamorph with periconiella-like branching patterns, and discusses the taxon's phylogenetic placement, and distinctive morphology. This discovery highlights the unexplored fungal diversity in East Asian forests.

The evolutionary transition to terrestrial life required overcoming several physiological hurdles; however, such challenges were amplified in desert environments. While several xeric-adapted arachnids utilize permanent burrows or "sit-and-wait" foraging strategies as possible energy conservation adaptations in harsh habitats, camel spiders exhibit a counterintuitive, high-energy lifestyle. To investigate the molecular underpinnings distinguishing Solifugae within Chelicerata, we utilized a comparative genomics framework that incorporates a newly sequenced, previously unpublished solifuge genome. We identified lineage-specific expanded orthogroups and evaluated selective pressures acting upon paralogous sequences within our ingroup solifuge species. Additionally, we also focused on fatty acid-associated proteins and heat shock proteins to elucidate how Solifugae may have evolved such anomalous behaviors compared to their arachnid relatives. Our analyses revealed significant signatures of positive selection within key gene families across the solifuge lineage. Notably, paralogs within the cytochrome P450 and biotinidase families showed consistent evidence of selection across all three taxa, suggesting specialized metabolic or detoxification requirements. Furthermore, we identified candidate loci implicated in axonal guidance and lipid metabolism, and a specialized fatty acid enzyme repertoire. While subsequent research is required to determine whether some of the genomic signatures unveiled here are shared across a broader phylogenetic distribution within Solifugae, we establish a critical baseline for future functional validation.

Remote island systems with small landmasses and reliable estimates of human occupancy are ideal model systems to disentangle the roles of global climatic changes and local human occupation on biota. Here, we used mitochondrial and nuclear genomic data from five endemic Fijian Lasioglossum (Homalictus) bee species to infer changes in effective population size (Ne). These ground-nesting bees are native, with non-specialised floral visitation habits, and distributed across the elevational gradient. All lowland species and populations showed strong signals of increasing Ne that correspond to the timing of human occupation of Fiji, but not Holocene climatic change. Highland populations, with greater isolation and present in regions less affected by anthropogenic impacts, did not show evidence of recent rapid increases in Ne. Population expansion rates across the elevational gradient differed between taxa, with significantly earlier and larger increases in predominantly lowland species than those with more restricted ranges in the highlands. This is consistent with the movement of people inland from coastal regions and into montane elevations of the island, and corresponding landscape changes that benefit the ecology of these bees. Specific life history traits of these bees, combined with substantive clearing of forest cover and floristic changes at lower elevations, has likely increased nesting opportunities and abundance of invasive floral resources. Our findings contrast with recent evidence that human occupation of Fiji has resulted in decreased ant biodiversity and raise the paradoxical possibility that human-mediated environmental changes may benefit some native montane tropical insect faunas.

Biological invasions on islands provide a natural framework to study how dispersal and connectivity influence evolutionary and ecological processes. The avian nest parasitic fly, Philornis downsi - first recorded in Darwin's finch nests in 1997 - causes high mortality in endemic land birds, yet its inter-island and sex-specific patterns of dispersal and genetic structure remain poorly understood. We use low-coverage whole genome sequencing to investigate genome-wide patterns of genetic diversity, directional migration and effective population size in P. downsi across five major Galapagos Islands and its native range in mainland Ecuador. We find evidence for a genetic bottleneck in the Galapagos, isolation by distance, and evidence that the island closest to the Ecuadorian mainland, San Cristobal, is genetically divergent from the other four islands sampled, despite retaining the highest genetic diversity. No evidence was found for sex-biased dispersal; however, sex-biased genetic structure was detected using only markers from inferred autosomal scaffolds. We found asymmetric gene flow with higher migration rates from San Cristobal westward to the other islands, matching the direction of both southeast trade winds and major cargo shipping routes. Our results suggest both natural and human-mediated colonisation of P. downsi from the mainland through San Cristobal to the other islands, followed by high inter-island dispersal among closely situated sink islands. Our findings are critical for prioritising islands for control strategies that will reduce P. downsi impacts on vulnerable endemic birds and underscore the value of understanding directional migration patterns for managing invasive species in metapopulations.

Morus rubra is native to the eastern United States, with its range extending into the Upper Midwest and southern Ontario, Canada. Its present distribution suggests that past glacial events in North America may have influenced the genetic structure of populations at the species northwestern range boundary. This study assessed genetic variation among six M. rubra populations believed to have experienced postglacial colonization using published nuclear microsatellite markers and sequences from the chloroplast trnL-trnF region. Five nuclear microsatellite markers previously developed for M. alba were successfully transferred to M. rubra, while the chloroplast trnL-trnF region provided an additional marker for evaluating chlorotype diversity. Nuclear microsatellite diversity was higher in southern unglaciated populations than in northern glaciated populations, a pattern consistent with the observed distribution of chlorotype diversity. Together, these results support ancient founder effects associated with leading-edge expansion following glacial recession and suggest that postglacial colonization contributed to the present-day genetic structure of M. rubra at its northwestern range boundary. Because M. rubra hybridizes with the naturalized invasive M. alba, reduced genetic diversity in marginal populations may increase their vulnerability to genetic swamping. The markers characterized in this study provide useful tools for population genetic research in Morus, and the findings have important implications for the conservation and management of marginal and threatened M. rubra populations in the Upper Midwest.

Ectomycorrhizal (EM) fungi are vital to forest ecosystems, supporting tree growth and survival. However, their inclusion in conservation policy and action remains limited and little is known about the status of their genetic diversity, which is essential for their long-term survival and adaptation. The Global Biodiversity Framework adopted a genetic indicator based on the effective population size, Ne, to monitor genetic diversity in all species. To date, it is still uncertain how Ne, a key parameter, can be reliably assessed in species with complex life history traits. Ectomycorrhizal fungi are a highly diverse group of taxa displaying haplodiplontic life cycles with partially clonal reproduction. Here, we review the literature to understand how these life history traits might affect Ne and its estimation in six species of EM fungi. We estimated Ne in 19 populations using eight genetic and genomic datasets from selected studies. We compared Ne estimates using Linkage Disequilibrium (LD) and Sibship Frequency (SF) methods. We tested how Ne estimates change due to partial clonality and genetic structure gradients and whether the number of genetic markers influence the precision of the estimates. We show a systematic bias in Ne estimations when large clones are present and when populations are not correctly delimited. We found both methods are not robust to these factors, which makes them unreliable for conservation assessment purposes in EM fungi. This study provides new perspectives for further research into the links between life history traits and the effective population size of ectomycorrhizal fungi.

Color polymorphism can facilitate local adaptation, maintain intraspecific diversity, or reflect early stages of speciation. Clownfishes (Amphiprion spp.) typically display a simple black, orange, and white pattern, but the saddleback clownfish (Amphiprion polymnus) shows striking variation in melanism and the number of vertical bars, which are thought to play a role in species recognition. In 2024, a revision on iNaturalist split A. polymnus into multiple species based solely on color pattern and geographic range. This raises the question of whether these morphs represent true species or intraspecific polymorphism, which we tested using genomic and image-based data. We sampled 97 individuals from seven populations across the species range and quantified color patterns from standardized photographs. Phenotypic and genomic analyses reveal a complex pattern of divergence. Image analysis identified three distinct phenotypic clusters, with A. polymnus, A. annamensis, and A. laticlavius each showing consistent differences in saddle shape and vertical bar extent. ADMIXTURE resolved three distinct genetic groups corresponding to the morphs. Pairwise FST (0.54-0.71) and dxy indicate extremely high differentiation between A. polymnus and A. annamensis, consistent with species-level divergence, whereas A. laticlavius shows much lower differentiation from A. polymnus (FST 0.09-0.18) and higher differentiation from A. annamensis (FST 0.64-0.66). Overall, phenotypic and genomic data show structured variation, but the status of A. laticlavius remains ambiguous. Our study reveals clear and structured divergence across the full range, yet the taxonomic interpretation of this variation remains inherently challenging. The key question remains: do these patterns reflect a single polymorphic species or a complex of closely related species?

Limb loss in vertebrate animals is surprisingly common despite imposing strong functional constraints. These pressures funnel species towards regions of limited ecological and phenotypic space. To date, snakes have been considered unique in having escaped this pattern. Using a new species-level phylogeny and comparative morphological and dietary datasets, we show that pygopods, a group of limbless Australo-Papuan geckos, have undergone a similar evolutionary trajectory to snakes. Our analyses provide evidence of exceptional morphological and diet evolution. This is exemplified by strong niche partitioning among genera through dietary specialization and greater than expected dietary disparity. Diversification in pygopods has also been driven by extreme phenotypic evolution, with pygopods encompassing much of the morphological space covered by all other limb-reduced lizards. Interestingly, the diversification of pygopods has resulted in only a modest number of species, emphasizing the decoupling of diversity and richness possible in adaptive radiations.

While asexual species can often outcompete their sexual counterparts over ecological timescales, their long-term evolutionary success is hindered by a diminished ability to purge deleterious mutations and to adapt to changing environments. However, some asexual species persist for millions of years, and a major question in evolutionary biology is how they do so. One solution is to occasionally reproduce sexually, as has been shown in a handful of primarily asexual species. Here, we investigate the possibility of rare sex in the clonal raider ant, Ooceraea biroi. We report the whole-genome sequence of a previously uncharacterized clonal line and, using population genetic and phylogenetic analyses, show that it originated through sexual reproduction between two extensively studied clonal lines. The mitochondrial genome of this clonal line differs from that of the maternal clonal line at only a single nucleotide, suggesting that the sexual reproduction event occurred within the past few hundred years. These results demonstrate that sex occurs sporadically in the clonal raider ant, allowing it to generate new genetic combinations and potentially to overcome some of the costs of asexuality.

Whole-genome alignment (WGA) is widely used for genome-scale phylogenetic inference, and most scalable WGA pipelines rely on progressive alignment guided by a pre-specified tree. Among progressive whole-genome aligners, Progressive Cactus is a successful state-of-the-art method. However, analyses of real and simulated avian data indicate that guide-tree choice can influence downstream tree inference; star guide trees do not remove this effect and can exacerbate long-branch attraction artefacts. We have developed a consensus strategy based on the Progressive Cactus framework by generating a small set of alternative guide-tree alignments and retaining only homology relationships consistently recovered across all alignments. In simulation experiments, consensus alignments improve precision, bring inferred site-pattern frequency distributions closer to those of the true alignments, and recover more true splits than single guide-tree alignments. In a real landbird (Telluraves) dataset, we observe a strong bias towards single binary guide trees and long-branch attraction for less resolved trees. While the reconstructed tree still depends on the phylogenetic method and taxa sampling, our consensus alignment has no clear bias. We implemented a hierarchical consensus workflow that only locally resolves uncertainty in the guide tree. Therefore, the computational cost increases only moderately, for example by an estimated 68 percent for a recently published large-scale alignment of more than 300 modern birds (Neoaves) taxa.